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Abstract 

Recent data of lattice measurements of the gauge-invariant nonlocal scalar quark con- 
densates are analyzed to extract the short-distance correlation length, 1/Ag, and to con- 
struct an admissible Ansatz for the condensate behaviour in a coordinate space. The cor- 
relation length values for both the quenched and full-QCD cases appear in a good agree- 
ment with the well-known QCD SR estimates of the mixed quark-gluon condensate, 2Ag = 
{q{igcy fj,uG pi_y)q) / {qq) = 0.8 — 1.1 GeV^. We test two different Ansatzes for a nonlocal quark 
condensate and trace their influence on the twist-2 pion distribution amplitude by means of 
QCD SRs. The main features of the pion distribution amplitude are confirmed by the CLEG 
experimental results. 
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1 Introduction 



The results of long-awaited lattice measurements of gauge-invariant nonlocal quark condensates 
have been published recently in These data provide a possibility to examine directly the 
models of nonlocal-condensate coordinate behaviour. Different models have been suggested in 
the framework of QCD sum rules (SR) for dynamic properties of light mesons 0, ^ ^, ^, |^. 
We used an extended QCD SR approach with nonlocal condensates || as a bridge to 

connect meson properties, form factors, and distribution amplitudes with the structure of QCD 
vacuum. For light meson phenomenology, the value of short-distance correlation length / in QCD 

vacuum, / ~ — , has a paramount significance, while the details of a particular condensate model 

are of next-to-leading importance. We demonstrate that the original lattice data in [|I| allow one 
to extract a reasonable value of the correlation scale being in agreement with our vacuum 
condensate models. At the very end, the lattice data support our conclusion on the shape of 
the pion distribution amplitude 0, ||, |^, and vice versa, the CLEO experimental data on 
pion photoproduction confirm our suggestion about the range of the correlation scale values in an 
independent way []TU| . 



2 Models of nonlocal quark condensates 

Small behavior. Let us start with the general properties of gauge-invariant nonlocal quark 
condensates (NLCs) Fsy{x'^) following from their definitions 



(:g.(0)E(0,a;)g,(x):) = M 



Fs{x') - '-^Fy{x'^ 



E{0,x) = Pexp -z / A^{y)dy'^ , (1) 



pa 



where a, p are spinor indices, the integral in the Fock-Schwinger string E{0,x) is taken along a 
straight-line path. First of all, the condensates Fsy^x"^), . . . should be analytic functions around 
the origin, and their derivatives at zero are related to condensates of the corresponding dimension. 
Expanding Fsyix"^) in the Taylor series at the origin in the fixed-point gauge A^{y)y'^ = (hence 
E{0,x) = 1), one can obtain 



(g(O)g(x)) _ (g(l + ^^Df + j,{xD)' + ■■■)?) 



(qq) {qq) 



6 / 9\ 2 



^ , 4(g(0)i-^(i-)) _ ^4(# ((xD) + ^,(xDf + ^(xDf + ■ ■ ■)q) 
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{qq)x^ {qq)x 

~ 3!3Q^ W 5!3g^ W + ^^"^ ^ " 

Here = d^j_ — igA^^t°- is a covariant derivative; rrig, the quark mass; J^, quark vector current; and 
the expansion coefficients appearing in @-(|]) are vacuum expectation values (VEVs) of quark- 
gluon operators Oi of dimension z, = (Oj). These expansions with the explicit expressions for 
the condensates have been derived in ||Tl|, (see Appendix A). Condensates of lowest dimensions 



Q^ = {qq), Q' = ig{q{G,,a,.)q)^ml-Q\ ^ , (4) 

form the basis of standard QCD SR [T^ and have been estimated, while the higher-dimensional 



VEVs, QJn),Ql^ )^Q%) ■ ■ ■ yet unknown. Note here that the matrix element of the 4-quark 
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condensate is not known independently. Instead, it is usually evaluated in the "factorization 
approximation", the accuracy is estimated to be about 20% [1T2| , 



{qViqqT2q) 



Tr (FiF, 
144 



■{qq? 



The "mixed" condensate is expressed in the chiral limit as rn^ ■ or 2Ag ■ 



\x=0 



4Q3 



rrin 



A 



2 ' 



(5) 



(6) 



A^ 

and the parameter fixes the width of Fs{x'^) around the origin. This important quantity has 
been estimated within the QCD sum-rule approach 



A 



0.4±0.1GeV2 |3 
0.5 ± 0.05 GeV^ 



(7) 



Estimates of A^ from instanton approaches |]T5l are somewhat larger: A^ > 2/p^ ^ 0.6 GeV^ where 
Pc ~ 1.7 — 2 GeV~^ is an average characteristic size of the instanton fiuctuation in the QCD 
vacuum. Taking into account all these estimates, in the following, we put a rather wide window 
0.6 > A2 > 0.4 GeV2 for its value ("QCD Range" in Figs.|, |). 

Large asymptotics from HQET. The large- |x| properties of the NLC Fsix"^) have been 
analyzed in detail in [|TB[ in the framework of QCD SRs for heavy-quark effective theory (HQET) 
of heavy-light mesons. It was demonstrated that for a large Euclidean x, NLC is dominated by 
the contribution of the lowest state of a heavy-light meson with energy A = (Mg — mQ)\mQ^oo, 
and Fslx"^) ~ exp{ — |a:;|A} (numerically, A is around 0.45 GeV). In the following, we shall take 
this asymptotic behavior for NLC at large |x| 



(8) 



Hints from QCD SRs, Gaussian Ansatz. To relate the NLC behaviour with the properties 
of mesons via QCD SR, it is convenient to parameterize the x^-dependence of (0)-(§|) by the 
distribution functions fx{c() o,'la a-representation for a propagator 



Fsy{x^)=J e-'^'"'/^ fs,v{a)da, where J fsy{a) da = 



1, S-CcLSG^ 

0, V-case, chiral limit. 



(9) 



Here we use the Euclidean interval = —x'^ > 0, and the subscript E will be omitted below 
for simplicity. The representation (^ allows one (i) to involve smoothly NLCs into diagrammatic 
techniques, and (ii) to clarify the physical properties of NLCs. Indeed, functions fs,v,..X'^) intro- 
duced in describe the distribution of quarks over virtuality a in nonperturbative vacuum. The 
moments of fsy,... coincide with Taylor expansion coefficients in (@)-(0). For example, we have in 
S case in the chiral limit 



fOO 



afs{a) da 



_ 1 [ {■■qP'q--) 
{■.qq:) 



POO 



a'^fs{a) da = 



(0) 



18Q3 



(10) 



with Ag meaning the average virtuality of vacuum quarks. Higher moments of /5(a) are connected 
with higher-dimensional VEVs. The difference between the Ansatzes for Fs{x'^) looks more pro- 
nounced just for its /5(a)-images. It is evident that distributions extremely concentrated at the 
origin fsia) ~ 5{a), S^^\a), . . . correspond to separate terms of the Taylor expansion (0)-(^. At 
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the same time, /5(a) - 
The Gaussian Ansatz 



CONST simulates a free propagation with zero mass: ~ CONSt/x^. 



2 

exp I ——X 



(11) 



takes account of a single scale - "inverse width" Ag of the coordinate distribution and corresponds 
to the virtuality distribution fs = S{a — Xg/2). It fixes only one main property of the nonper- 
turbative vacuum — quarks can flow through the vacuum with a nonzero momentum k, and the 
average virtuality of vacuum quarks is (/c^) = \q/2. The Gaussian behavior at very large does 
not correspond to the expected NLC asymptotics (H). But for the moment QCD SRs, that deal 
with the smearing quantities - moments of distribution amplitudes p, ^ , form factors |], |1^ , 
the incorrect asymptotics of NLC as well as subtle details of the Ansatz shape are expected to be 
not so important. It is interesting that the Gaussian behaviour is supported by a model of the 



nonperturbative propagator, [T^, based on simple "local-duality" arguments. Namely, Eq.(16) in 
jT^I demonstrates a behaviour rather close to Eq.([TT|) within a physically important region 1 Fm. 
Unfortunately, corresponding Fs{x'^) decays too quickly beyond this region, it becomes negative 
and oscillating, which is physically unclear. 

Certainly, a more realistic model of /(a) should possess a finite width: we expect that it is a 
continuous function concentrated around a certain value Ag/2 and rapidly decaying to zero as a 
goes to or 00. Moreover, the continuous distribution /5(a) over virtuality is directly related to 
the pion distribution amplitude ipn{x) (for details, see sect. 4), as it was demonstrated with the 
help of the "nondiagonal" correlator in (^,0. 



\5{a) \S{a- Xl/2) 




Figure 1: Graphical illustration of different NLC models in a-representation. Thick bell- like (blue in a colour 
view) lines represents the advanced Ansatz (solid line) in comparison with fs{a;5) (dashed line), vertical line 
(green in a colour view) at a = - Gauss Ansatz, and vertical line (red in a colour view) at a = - the local 
limit to standard condensates. For comparison, we depict here the free propagation a-distribution in a dashed thin 
line. 

"Advanced Ansatz" . To construct models of nonlocal condensates, one should satisfy some 
constraints. For instance, if we assume vacuum matrix elements {:q{D'^)"^q:) to exist, then the 
function f{a) should decay faster than l/a"^~^^ as a ^ 00. If for all m, all such matrix elements 
exist, a possible choice could be a function /(a) ~ a" exp(— aa), etc. at large a. The opposite, 
small-a, limit of f{a) is determined by the large-|x| asymptotics (P) of the function Fs{x'^). This 
means that /(«) ~ exp(— A^/a) in the small-a region. By the simplest composition of both the 
asymptotics |^ 0, we arrive at the class of Ansatzes 0, |^ 

that gives a coordinate behavior 

F ix- n) - ^" (^^^^^) ( 2. y 
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where Kn{z) is the modified Bessel function. The distribution fs^cx) = fs{(y;n = 1) with the 
iq ~ 0.45 GeV and cr^ 



parameters ~ 0.45 GeV and cr^ ~ 10 GeV ^ is presented in Fig.l in comparison with the 



distribution fs{(y',n = 5). In this case, the behavior of Fg{x^) = Fs{x^; 1) is similar to that of 
the massive scalar propagator with a shifted argument: 

_ i^i(Av/4o-2 + |x^|) 2a A K, (2Aa) 

K,{2Aa) ^4^2 + 1x2 1' ~ ^ a ' K, {2Aa) ' ^^^^ 

The short-distance correlation scale in (|^) appears to be equal to at Aa ^ 1, reproducing 
the single instanton-like result where the parameter a imitates the instanton size Pc- In the 
opposite case at Aa ^ 1, is proportional to 2A/a. The asymptotics of F^(x^) at large |x| ^ 2a 
I is determined by the parameter A, 

F#(x^)U^^oo-^W-^^e-^^ (15) 



where A plays the role of an effective mass. This interpretation is transparent in the momentum 
representation of the NLC: Fg{p^) tends to the usual propagator form at small values of all the 
arguments 

[a^W+W\) A A' 



2\ 

r\-/ 



' " {2Aa) ^A2 + |p2| ■ A2 + p2 

"Advanced" Ansatz has been successfully applied in the nondiagonal QCD SR approach to the 
pion and its radial excitations [0, and the main features of the pion have been described: the 
mass spectrum of pion radial excitations vr' and tt" in close agreement with experimental data. 



3 Vacuum correlation length from lattice data 

Here we consider the results of fitting the lattice simulation data for both the kinds of Ansatz 
F^'^ix?) introduced previously for the nonperturbative part of a correlator. 

Processing the Pisa lattice data. The correlator Cq{x) = - E/=i(Tr [g^(0)^(0, x)g-^(x)]) 
was measured on an Euclidean lattice fTJ at four points, x = 2a, 4a, 6a, 8a, inside 1 Fm, where 
a ~ 0.1 Fm is the lattice spacing. The simulation was performed with four fiavors [Nj = 4) of 
staggered fermions with the same mass rriq for every fiavor and SU{3) Wilson action for the pure 
gauge sector. Both the full QCD case (f-QCD), i.e. , including the effects of "loop" fermions with 
the flavour masses m^a = 0.01, 0.02 and the quenched case (q-QCD where the effects coming 
from loops of dynamical fermions are neglected) with a set of flavour masses m^a = 0.01, 0.02, 
0.05, 0.10 were considered, see Appendix B and [|| for details. 

The correlator can be written as a sum of a perturbative-like term, B ■ F|*(x^), proportional to 
niq/x'^ at very short distances, and a nonperturbative part, A ■ F^''*(x^). Parameters of the latter 
are the main goal of the flt. For the perturbative part, the authors of have used just the short- 
distance asymptotics; whereas for the nonperturbative one, an exponential Ansatz F^^^(x'^) = 
exp(— /io|a;|) that corresponds to the asymptotics of very large distances x'^ ^ 1 Fm^. In other 
words, these two approximations, involved in the flt simultaneously, are adapted for different 
regions of x. It seems more dangerously that the exponential Ansatz is a non-analytic function of 
x'^ at the origin: it depends on |x| = y/x^ and its derivatives with respect to x'^ do not exist at the 
origin. Therefore one loses any connection with local VEVs appearing in the OPE and with the 
corresponding interpretation of the correlation length, see section 2. Nevertheless, the flt of the 

^The asymptotics starts with 3> 2 Fm^ at the mentioned value of cr-parameter. 
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lattice data has been performed and demonstrated very small y? . Note, the values of l^d.oj. 
in the Tables below should be considered as purely indicative of the best fit quality. One could 
not interpret these values in a standard statistical sense, because their true norm is unknown, see 
brief discussion in |Q. The extracted quantities, A, /xq; are presented in Table 1 (for every 
run) for a comparison with our results. 



Table 1: Modifiecl|| fit QTp , Fs^\x) = exp (-|a;|/xo) 



/3, theory 


a ■ niq 


aMo X 102 




B 




5.35, f 


0.01 


0.49(13) 


0.16(4) 


1.3(3) 


1.3-10-2 


5.35, f 


0.02 


1.7(5) 


0.26(4) 


0.95(5) 


6.4-10"^ 


6.00, q 


0.01 


1.6(5) 


0.16(4) 


0.9(1) 


7.6 - 10-2 


5.91, q 


0.02 


2.3(7) 


0.26(3) 


1.25(7) 


5.2 - 10-2 


6.00, q 


0.05 


1.8(4) 


0.34(2) 


1.4(2) 


0.2 


6.00, q 


0.10 


5.6(5) 


0.55(1) 


1.0(1) 


1.3 - 10-2 



We test different NLC Ansatzes and extract the parameters of the correlator using three-step 
procedure. At the first step, we fit rearranged formulas for the correlator. Note that the masses 
of light flavors appear unnaturally large in the lattice simulation and cannot be considered small 
(for the values of masses in MeV, have a look at the axes of Fig. 3). For this reason, one should 
take account of all possible mass terms in both parts of Cq: 

C,{x) = B - Ff{x^) + A - (a;2), (16) 
B ~ Ff{x') = m, (^^K,{m,x)) ; F^ix') = F^^\x') + + .... (17) 

Following [Q, we fix the perturbative-like part in the form of free propagation F^* = —Ki{mgx) 
with mass nig. For nonperturbative part, we keep all the known mass-terms from the expansion 
(0), see the expression for Q^^ in Appendix A. 



Table 2: Fit for the Gaussian Ansatz, F^{x) = exp 



P, theory 


a ■ ITLq 


O^Aq X 102 


(a ■ Xl? 


B 


xVNd.o.f. 


5.35, f 


0.01 


0.28(4) 


0.10(3) 


1.67(14) 


0.075 


5.35, f 


0.02 


0.69(12) 


0.18(3) 


1.67(21) 


0.125 


6.00, q 


0.01 


0.93(14) 


0.10(3) 


2.19(65) 


0.01 


5.91, q 


0.02 


0.92(14) 


0.16(3) 


2.31(27) 


0.65 


6.00, q 


0.05 


0.60(7) 


0.23(2) 


1.7(6) 


1.78 


6.00, q 


0.10 








2.7 



Table 3: Fit for the "Advanced" Ansatz, F/(x),A = 0.5 GeV 



/?, theory 


a ■ niq 


a^Ao X 102 


(a - Xl? 


B 


xVN,,o.f. 


5.35, f 


0.01 


0.29(5) 


0.11(3) 


1.63(19) 


0.03 


5.35, f 


0.02 


0.76(16) 


0.23(5) 


1.57(30) 


0.05 


6.00, q 


0.01 


0.97(17) 


0.12(3) 


1.95(63) 


5 - 10-4 


5.91, q 


0.02 


1.04(20) 


0.21(5) 


2.12(40) 


0.3 


6.00, q 


0.05 


0.79(14) 


0.38(6) 


1.59(9) 


0.6 


6.00, q 


0.10 


2.72(32) 


1.41(7) 


1.13(6) 


0.3 



^Note here that the authors of used the dimensional coefficient Bq related to our dimensionless parameter, 
Bo = rrigB. 

"in the original table jl], the dimensionless combination aB^ has been shown. Due to the previous footnote we 
recalculate the corresponding B parameter for their fits. 
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As a result of the fit, we extract, from lattice data, an intermediate correlation scale Xl that 
parameterizes F^'^i^x"^) in ([T7|) instead of Xg, 



and depends on lattice conditions, namely m^, a, and /5. We expect that this quantity coincides 
in the chiral limit with A^, determined in (P) within the massless QCD, i.e. , X\{mq) Im^^o — * 
Aq = Ag(/i|) on the lattice normalization scale The results for A, B and A|(mg), obtained in 
various casesQ are collected in Tables 2 and 3; let us outline the main features: 

(i) The for the cases of both Ansatzes look higher than in Table 1, but still sufficiently 
small, especially for f-lattice. An exception is the last run for the q-lattice with the largest mass 
irtq ■ a = 0.10 that corresponds to ^ 200 MeV [|I|. To process the data with so a huge mass, 
one should involve a lot of mass-terms into the fit formula For this reason, we exclude this 
run results from subsequent analysis. 

(ii) The extracted coefficient B that fixes the perturbative-like contribution should not change 
significantly from one run to another for both Ansatzes and for both kinds of lattices. This 
property can signal about a good quality of the fits, it confirms the reliability of the fit at least 
for the f-lattice case, see Tables 1 and 2. 

(iii) The fitted parameter X\{mq) that fixes the behavior of the nonperturbative part has a 
strong and monotonic dependence on m^, in contrast with the parameter B, see item (ii). 

Extrapolation to the chiral limit. At the second step, we extrapolate the intermediate 
X\{mq) to the chiral limit, suggesting a simple linear dependence on 0. The linear law 

seems to be rather naive, but the observed strong dependence of X\{mq) on the quark mass is well 
supported by the data, see graphics in Fig.^, |^. Really, the linear extrapolation of the first three 
run results for the q-lattice is self-consistent for both Ansatzes (Fig.|^(b), ||(b)), and the results 
are in a reasonable agreement with the corresponding f-lattice results (Fig.0(a), ^(a)). 




10 20 30 40 50 20 40 60 80 100 



Figure 2: Graphical illustration of the fit results for \\ - black points with error bars; for the chiral limit ( 
mq 0) and the evolution to the scale 1 GeV^ for the Gaussian Ansatz both in full (part (a)) and quenched 
(part (b)) LQCD. Dashed black lines on both parts correspond to chiral limit procedure with taking into account 
error-bars, black blobs - to the central points of the chiral limit, and red blobs - to the resulting values of on 
the QCD scale of 1 GeV^. Red arrows show the overall error-bars of the extracted A^(/i^ = 1 GeV^), whereas green 
thick hnes bound the QCD preferred region. 

The explicit results for Aq are presented in the first line of Table 4: note, Aq for the q-lattice 
is not smaller than for the f-one for both Ansatzes; the large error bars for Ag appear due to the 
roughness of the chiral extrapolation procedure. 

**The Gaussian Ansatz has been tested in a lattice in j2l] (after our suggestion in a private communication) 
without any mass corrections. The fit-result, X\[mq ■ a = 0.01) — 0.46(5), appears somewhat larger than our 
fit-result in Fig.^(a). 
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Figure 3: Graphical illustration to three-step procedure for the case of Advanced Ansatz both in full (part (a)) 
and quenched (part (b)) LQCD. 



To clarify the reliability of the approximation, we repeat the same extrapolation procedure with 



the dimensionless "lattice quantities", L 



Xl) and M = a ■ niq. Lattice spacing a involved 



into this extrapolation is different for different runs (see Appendix B, Eqs.( [B.l]) -( [B.2| )). To return 
to "physical" quantity Aq at the very end, we adapt an average spacing corresponding to the 
q-lattice case, and a similar average one a/ - to the f-lattice case (see Appendix B). It is clear, 
that this procedure is even more crude than the previous one. Nevertheless, the corresponding 
Ao well agree with the previous "three-point" q-lattice result, compare the first and third lines in 
Table 4. But the procedure falls down for f-lattice data, the results appear too small. 

At the third step, we perform an evolution to continuum normalization scale. To 
compare the results for the shape of NLC ^^(x^,/!^) (or distribution /(a,/i^) at scale /i^) with 
the same quantity taken on another scale Q^, the corresponding evolution law in /i^ is required. 
For a general case it looks as ambitious and rather complicated problem that is yet unsolved. But 
in our partial case if we fix how A^(/x^) evolves with /x^, then the evolution law for both Ansatze 
is also fixed. Therefore, we shall consider the evolution of single characteristic of the shape - 
Ag(yU^), setting the equivalence A^(/i|) = Aq in a lattice. We live not in a lattice, but, instead, in 
a continuum. Therefore we need a procedure to relate lattice measurements with our observables 
in continuum, the corresponding one-loop evolution law for this kind of a lattice presented, e.g. , 
in 



22 



A;(/^') ■ [ln(/i/A 



MS, 



7a 



X2 
-^0 



■(3 + ln(AL/AAf5) + 6.54 



7a 
bo 



;i8) 



The continuum quantities enter into the l.h.s. of Eq. ([Iq) ; while the lattice quantities, into the 
r.h.s.. Here Kms/^l = 76.44, = 11/3Ca — 4/3TrA^j - first coefficient of the /9-function, and 
7a = (2C^ — SCp) = —14/3 - one-loop anomalous dimension for A^(yU^), calculated in [^. So, 
using formula (0) to return to continuum Ag(yU^ = 1 GeV^), we obtain the final results for A^ 
appearing just in the "QCD range", see Figs.^, ^ Let us consider these evolved results, presented 
in Table 4, in more detail. 

(i) The Gaussian Ansatz. The mean values of A^ for f- and q-lattices become closer to one 
another after the evolution starting from the corresponding different values of Ag. This "focusing" 
effect is due to a difference of the evolution in both the cases: for the f-lattice (/5 = 5.35; Bq = 
25/3, Nf = 4) the transition factor from Eq.(|18D is 2.88, while for q-lattice {p ^6; bo = 11, Nf = 
) - 2.11. Note here that if one attempts to exchange, by hand, these evolution laws (q for f and 
vice versa), then the final numbers diverge out of the QCD range. The observed focusing may 
demonstrate a complementarity of the chiral limit results to the evolution law. 

(ii) The Advanced Ansatz. The value of A^ for f-lattice, 0.5 GeV^, appears to be very close 
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to an average for the Gaussian Ansatz in the left side of Table 4. The result for the q-lattice 
is less than this average estimate and located near the low boundary of the QCD range. But, in 
virtue of huge error bars, the latter result does not contradict the f-lattice result. 

Finally we can conclude that Pisa lattice data reported in [l| really "feel" the short distance 
correlation scale in the QCD vacuum. The data processing explicitly demonstrates that the 
extracted mean values of are in agreement with the estimates from the QCD SR approach (|^) 
for all the considered cases. Moreover, the results are in agreement with the old lattice result, 
Ag ~ 0.55 GeV^ , obtained in p4| on a q-lattice. Huge errors bar are the problem of these 
estimates, and this does not allow one to confirm the agreement with the "QCD range" once and 
for all. The main uncertainty follows from the chiral limit procedure ("second step"). To reduce 
the uncertainty, one should improve the theoretical part in (|l^) as well as the "lattice" part of 
the fit. Namely, we need more numbers of the lattice runs with a moderate/small quark masses 
M = a - TUq for a reliable extrapolation; to process all existing results, we should include the most 
important subset of mass-terms into the nonperturbative part of the fitted formula (JT 

Table 4: On the analysis of different approaches to the chiral limit 





Full QCD 


Quenched QCD 


Full QCD 


Quenched QCD 


A^ [GeV^], Gauss- Ansatz 


Aq [GeV^], Advanced- Ansatz 


Dimensional 
Chiral Limit (Aq) 


0.19 ±0.10 


0.23 ±0.16 


0.175 ±0.1 


0.175 ±0.26 


Evolution to 
111 = 1 GeV^ 


0.55 ±0.28 


0.48 ±0.33 


0.50 ±0.27 


0.37 ±0.54 


Dimensionless 
Chiral Limit (Ag) 


X 


0.24 ±0.15 


X 


0.19 ±0.23 


Evolution to 
III = 1 GeV2 


X 


0.50 ±0.31 


X 


0.40 ±0.50 



Results of extrapolation to uiq —>■ are evolved to the standard normalization scale /Xq = 
1 GeV^. Crosses in full QCD columns mean the breakdown of chiral limit procedure. 

Another problem is to extract the quark condensate value \{qq) \ from the lattice QCD data. 
In the fit it is just the coefficient A divided by the number of flavors Nf = 4. But the three-step 
procedure fails for this datap] providing too small values for the quark condensate. 

Another possibility is to construct the RG-invariant quantity rnq{fj,'^){qq){fi'^) in the lattice to 
avoid both the chiral limit and the renormalization effects. In this way, we obtained different 
estimates for every run in full the lattice QCD; the estimate region is 



\2mq{qq) \ = (2 - 5) ■ 10"^ GeV^ , 
that should be compared with the well-known value fixed by the current algebra 



(19) 



\{m^uu) + {mjd)\ = -/>^(l±0(m^)) 



1.7-10-^ GeV^ 



(20) 



^^^^The authors of |Q] also did not obtain reasonable estimates for (qq) using this kind of data, therefore they have 
performed an individual measurement of the quantity, see Eq. (3.3-3.5) in 
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The estimate (|T9D produces for real QCD case with the current masses rriu^d — 5.5 MeV 



GeV- 



■2M1/3 



265 - 358 MeV vs the standard value 250 MeV, 



(21) 



The values are overestimated, although the lowest one corresponding to the run with niq ■ a = 0.01 
looks reasonable. 

4 Nonlocal quark condensates and pion distribution 
amplitude 

The pion distribution amplitude (DA) of twist-2, y9^(a;, /x^), is a gauge- and process-independent 
characteristic of the pion that universally specifies the longitudinal momentum xP distribution of 
valence quarks in the pion with momentum P (see, e.g., for a review). 



(0 I d{0)-f>'-f5E{0,z)u{z) I 7r(P)) 



=0 



Jo 



V?,r(a;,/i 



(22) 



Due to factorization theorems [26, 27|, it enters as the central input of various QCD calculations 
of hard exclusive processes. Here we illustrate how a value of the correlation scale (~ 1//^) can 
affect the shape of the pion DA. First, we consider NLC QCD SR for DA moments that provides 



the smearing quantities, moments 



{2x — 1) (pTr{x)dx, to restore a profile (fnix) of the 



pion DA. The NLC QCD SR are based on different kinds of Gaussian Ansatzes p|, ^ for NLCs 
that naturally appear in the theoretical (r.h.s.) part of the SR. 

Pion DA from NLC QCD SR for pion DA moments. The Gaussian Ansatz. The SR 

involves 5 different kinds of nonlocal condensates in addition to the scalar condensate contribution, 
A$s(a;; M^), for details see [El 11, H. The scalar NLC contribution results from the "factorization 



2.5 



1.5 



0.5 



(^i{x) at A2 = 0.4 GeV2 



(a) 




2 . 5 



1.5 



. 5 




Figure 4: Graphical representation of (p^i(x, fi^) (part (a)) and 93^2(2;, A^^) (part (b)) at the characteristic renor- 
malization point /i^ ~ 1 GeV^. The thick sohd hnes in both plots denote iy9°^*2(^)i i-^-; the best fit to the determined 
values of the moments, whereas dashed lines illustrate admissible options approximately with < 1. 



approximation" to the nonlocal four-quarks condensate, and its accuracy is yet unknown, compare 
with the approximation, Eq.(|^). But in any case the contribution A$5'(x; M^) is numerically the 
largest one for not too high moments. Therefore, the main features of the shape of ^t^{x) in Fig.^ 
is, roughly speaking, the net result of the interplay between the perturbative contribution and the 
nonperturbative term A$5(x; M^) that dominates the r.h.s. of the SR. 

The QCD SR predicts the values of moments {C,^)n{f^'^ — 1 GeV^) within their error bars. For 
this reason, one obtains, after restoration a "bunch" of admissible DA profiles corresponding 
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to the moment error bars, rather than a single sample of profile. These profiles are shown in Fig.^ 
by dashed lines, in addition to the optimal one (thick solid line) that corresponds to the best fit at 
xlpt ~ 10~^. Comparing DA profiles at different values of in Fig.^(a) and (b), one can conclude 
that the larger is the correlation scale the smaller is the concavity in the middle of the profiles, 
and their shape becomes closer to the shape of "asymptotic" DA ipas{x) = 6x(l — x). Therefore, 
a trial bunch (is not shown here) corresponding to the value = 0.6 GeV^ at the boundary of 
the introduced QCD range contains mainly convex in the mid-point profiles that are close to the 
asymptotic one. 

We have established in |I0| that a two-parameter model <^7r(2^; ^2, 0.4), the parameters being 
the Gegenbauer coefficients 02 and 04 (as also used in [^]), enable one to fit all the moment 



constraints for {^'^)n- For completeness we write explicit formulae for the optimal DA models 

<^°Pt(a;) = ^^{x)\l + aT^'^ ■Cl^\2x-l)+aT^^^ ■Cl^\2x-1)\ , (23) 
a°P*^(/i2) = +0.188 , aT\fi^) = -0.130 , at Aj = 0.4 GeV^; 
a°P*2(^2^) _ +0.126 , al^'^fi^) = -0.090 , at Aj = 0.5 GeV^ 

at /i^ ~ 1 GeV^. In this way, the admissible bunches of profiles can be mapped into the 
a2(/i^), a^ifi"^) plot and then can be evolved to a new normalization point [|1^, [iskYi see 
slanted rectangles in Fig.^. 

Pion DA vs CLEO data. Recently, the CLEO collaboration measured the 7*7 tt" 
form factor with a high accuracy. These data sets were processed by Schmedding and Yakovlev 
(S&Y) using a NLO light-cone QCD SR analysis. They obtained the constraints to (02, 04) 




Figure 5: The parameter space of (02, 04) pairs, corresponding to the allowed values of the second and fourth 
Gegenbauer coefficients: (i) calculated within the NLC-SR approach for three different values of and evolved 
to ^ig^Y — 5-6 GeV^; (ii) confidence regions extracted by Schmedding and Yakovlev Q from CLEO data 
Contour lines show 68% (solid line) and 95% (dashed line) confidential regions, ft is shown how our estimate of 
the confidence region with = 0.4, 0.5, 0.6 GeV^ overlaps with those displayed in Fig. 6 of Ref . . Bold dot in 
the plot marks the parameter pairs for the asymptotic DA, full square - for Chernyak-Zhitnitsky DA. 

coefficients in the form of 95% (2o"-deviation criterium) and 68% (la-deviation criterium) confi- 
dence regions, see ellipse-like contours in Fig. ^ with the central point = 0.190, = —0.140 
that (approximately) corresponds to an average normalization point of /x^^y = 2.4 GeV. 

The regions enclosed by the slanted rectangles bounded by the short-dashed line and solid 
line in Fig.^ correspond to the bunches of DA displayed in Fig.^(a) and (b). The central point 
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02 = 0.142, 04 = —0.087 corresponding to the optimum profile of tlie v?i(a;)-buncli (A^ = 0.4 GeV^ 
Fig. 3 (a)) at definitively belongs to the central S&Y region. The (/92(x)-bunch is, however, 

mostly outside the central 68% region though still within the 95% confidence region. Finally, the 
third slanted rectangle limited by the long-dashed contour and shifted to the upper left corner of 
the figure in Fig.|^ corresponds to the trial bunch of NLC-DA with = 0.6 GeV^. This value 
falls actually outside the standard QCD NLC-SR bounds in Eq.(^ for A^. Remarkably, the image 
of this region in Fig.^ lies completely outside the central region as a whole. Therefore, we may 
conclude that the CLEO data prefer the values A^ = 0.4, 0.5 GeV^ and probably do not prefer 
the value A^ = 0.6 GeV^, in full agreement with previous QCD SR estimates. Now the conclusion 
is supported by the lattice results presented in section 3. The quantitative details of the above 
qualitative discussion can be found in [|10| . 

Pion DA from nondiagonal correlator. The "Advanced" Ansatz. An approach to 
obtain directly the forms of the pion and its first resonance DAs, by using the available smooth 
Ansatz for the correlation functions /(a) of the nonlocal condensates, was suggested in papers 

01 . The sum rule for these DAs based on the nondiagonal correlator of the axial and 

pseudoscalar currents has the vanishing perturbative density. 



+ . . . 



M2 



1 — X 



X2 



2M2) 



fs{xM') +\x + 



X2 



2M2) 



fs{{l-x)M^ 



(24) 



only quark and mixed condensates appear in the theoretical part of the SR 0. The SR results in 
the elegant Eq. ( p4|) from the approximations both in the theoretical (for a detailed discussion, see 
0) and the phenomenological parts (see 0). In virtue of the approximations, a single correlation 
function /5(a) appears in the r.h.s. of Eg . (plf) that determines the profile of '^t^{x). 

Equation ( p4D demonstrates, in the most explicit manner, an important relation: the distribu- 
tion ^Pt^{x) of quarks inside the pion over the longitudinal momentum fraction x (on the l.h.s.) is 
directly related to the distribution fs{ot) over the virtuality a of the vacuum quarks on the r.h.s.. 
Note that a similar relation was obtained in an instanton-induced model |2^ . 

The approaches to extract ^-ni^x) with the help of SR (|^) were discussed in detail in 0, |^, 
here we concentrate on the final result for the profile in the case of the Advanced ansatz ([1^ . 
Note only that these approaches do not provide the behaviour of the profile in the vicinity of the 
end points, the reliable predictions are expected around the mid-region. 




Figure 6: Graphical representation of ipTr{x) in the 3-resonance approximation. Sohd hne corresponds to AP 
0.4 GeV^, short-dashed hne - A/^ — 0.6 GeV^, long-dashed line - Af^ = 1 GeV^, fine line - asymptotic DA. 



The shape of ^T^i^) strongly depends on the spectrum of the resonances in the l.h.s. of Eq.(^) and 
weaker - on the value of A^. For the model of equidistant, infinite number of narrow excitations 
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like the "Dirac comb" we have obtained the profile, very close to the asymptotic DA 0]. If we 
choose the spectral density in accordance with the current knowledge^, i.e. , containing only 3 
radial vr-excitations with the masses m^, ^ 1.3^ GeV^ , m^„ k. 1.8^ GeV^ , m\,„ ~ 4.7 GeV^ , 
then we obtain the set of admissible profiles presented in Fig.^. It is naturally to suggest that 
^Tr{x) is situated between these two possibilities. This result qualitatively agrees with the profile 
behaviour for the bunches ipi{x) and (x) in Fig.^. 



5 Conclusion 

We consider the admissible Ansatzes for the coordinate behaviour of the scalar quark nonlocal 
condensate being of importance in QCD SRs. These Ansatzes depend on the parameter — the 
short-distance correlation scale in the QCD vacuum that controls the corresponding coordinate 
behaviour. We analyse the lattice simulation data for the scalar NLC from and test two 
different Ansatzes. The correlation scales are extracted following the procedure that includes a 
fit of lattice data, an extrapolation to the chiral limit, and an evolution of the obtained lattice 
results to the characteristic scale /i^ = 1 GeV^ of a continuum QCD. 

The scale thus extracted does not visibly depend on the kind of a lattice (full or quenched) 
QCD, as well as on the kind of the Ansatz. The value of the scale appears in a good agreement 
with the QCD range A2(1 GeV^) = (0.4 - 0.55) GeV^, see Figs.2-3 and Table 4 in section 3. The 
agreement looks unexpectedly good in view of a roughness of the mentioned procedure. 

The scalar condensate and the scale A^ substantially determine the shape of the pion distri- 
bution amplitude by means of QCD SRs. Both kinds of the considered Ansatzes lead to similar 
shapes of the pion DA. The pion DA following from QCD SR for the moments p|, |l^] is rather 



sensitive to the value of A^. The bunches of DA profiles corresponding to the A^ values from the 



QCD range do just agree with the constraints that follow from the CLEO experimental data 

This and previous @, |^, |^, |^, |10| considerations demonstrate a close link between the correlation 
scale in the QCD vacuum and the shape of the pion DA. 
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Appendix 
A Basic set of local condensates 



Here we write down the explicit expressions ||TT| for the condensates that determine the first 
and second Taylor expansion coefficients in (El), 



''Only tt' and tt" are well established, while the mass of vr'" is fixed in analogy with other light meson spectra. 
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Q(o) = ^Ql-i^Ql-^Ql + Ql (A.i) 

Q\^) = + 3mQ5 _ 3^3g3^ (A 2) 

Qfo) ~ ~ 1^2 + ^Qs ~ 1^4 + SQg, (A. 3) 

g[^) = 5(3gj-g2'-3g5 + gD-i5m(Q' + ^Q'-"^'Q'), (a.4) 

where the quark condensate basis was chosen in the form 

Q' = {qq), Q' = ^9{q{G,.a,M, Q' = {t^J'/4), (A.5) 
for the condensates of lowest dimensions, and 

Ql = <qG^^G^^q>, Ql = i < qG^^G^^'jr,q>, (A.6) 
Ql = < qG^xGx^a^yq >, Ql = i < qD^J^a^^q >, 



for the condensates of dimensions 7. The basis elements of dimension 8, and A, enter into the 
expansion only of the vector NLC that is not analyzed here. For this reason, we do not show 
it here and refer the reader to article [|ri[], Eqs.(3.10-3.11). 

B Details of Pisa lattice simulations 

For the full QCD case (with dynamical fermions), the nonlocal condensates were measured on a 
16^ X 24 lattice at = 5.35 (/3 = Q/ g^, where g is the coupling constant) and two different values 
of the dynamic quark mass: a ■ rriq = 0.01 and a ■ rriq = 0.02 with the following values for the 
lattice spacing: 

a{p = 5.35) ^ 0.101 fm, for a ■ rUg = 0.01 ; 

a(/? = 5.35) ~ 0.120 fm, for a ■ = 0.02 . (B.l) 

For the quenched QCD case the measurement was performed on a 16^ lattice at (3 = 6.00, with 
the quark mass a ■ ruq = 0.01 for constructing the external-field quark propagator, and also at 
/3 = 5.91, with the quark mass a ■ rriq = 0.02. In both the cases, the value of (3 was chosen in order 
to have the same physical scale as in full QCD at the corresponding quark masses, thus allowing 
a direct comparison between the quenched and the full theory. In the quenched case, the lattice 
spacing is approximately 0]: 

a^^^^\p = 6.00) ~ 0.103 fm, for a ■ = 0.01, 0.05, 0.1 ; 

a(^*^)(/3 = 5.91) ~ 0.120 fm, for a ■ = 0.02 . (B.2) 

References 

[1] M. D'Elia, A. Di Giacomo, and E. Meggiolaro, Phys. Rev. D59 (1999) 054503. 

[2] S. V. Mikhailov and A. V. Radyushkin, JETP Lett. 43 (1986) 712; 
Sov. J. Nucl. Phys. 49 (1989) 494. 

[3] A. P. Bakulev and A. V. Radyushkin, Phys. Lett. B271 (1991) 223. 

14 



[4] S. V. Mikhailov and A. V. Radyushkin, Phys. Rev. D45 (1992) 1754. 
[5] A. V. Radyushkin, |arXiv:hep-ph/9406237| . 



[6] A. R Bakulev and S. V. Mikhailov, JETP Lett. 60 (1994) 150; 
in *Vladimir 1994, Quarks '94* 574-580. 

[7] A. P. Bakulev and S. V. Mikhailov, Z. Phys. C68 (1995) 451; 
Mod. Phys. Lett. All (1996) 1611. 

[8] A. P. Bakulev and S. V. Mikhailov, Phys. Lett. B436 (1998) 351. 

[9] J. Gronberg et al, Phys. Rev. D57 (1998) 33. 
[10] A. P. Bakulev, S. V. Mikhailov, and N. G. Stefanis, Phys. Lett. B508 (2001) 279. 
[11] A. G. Grozin, Int. J. Mod. Phys. AlO (1995) 3497. 

[12] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl. Phys. B147 (1979) 385; 448; 519. 

[13] V. M. Belyaev and B. L. loffe, Sov. Phys. JETP 57 (1983) 716; 

A. A. Ovchinnikov and A. A. Pivovarov, Sov. J. Nucl. Phys. 48 (1988) 721. 

[14] A. A. Pivovarov, Bull. Lebedev Phys. Inst. 5 (1991) 1. 

[15] A. E. Dorokhov, S. V. Esaibegian, and S. V. Mikhailov, Phys. Rev. D56 (1997) 4062. 
M. V. Polyakov and C. Weiss, Phys. Lett. B387 (1996) 841. 

[16] A. V. Radyushkin, Phys. Lett. B271 (1991) 218. 

[17] S. V. Mikhailov and A. V. Radyushkin, Sov. J. Nucl. Phys. 52 (1990) 697. 
[18] A. G. Grozin and Oleg I. Yakovlev, Phys. Lett. B291 (1992) 441. 
[19] A. A. Pivovarov, Phys. Atom. Nucl. 59 (1996) 891. 

[20] H. G. Dosch, M. Eidemiiller, M. Jamin, and E. Meggiolaro, JHEP 07 (2000) 023. 

[21] E. Meggiolaro, Nucl. Phys. Proc. Suppl. 83 (2000) 512. 

[22] R. Altmeyer et al, Nucl. Phys. B389 (1993) 445. 

[23] A. Yu. Morozov, Sov. J. Nucl. Phys. 40 (1984) 505. 

[24] M. Kremer and G. Schierholz, Phys. Lett. B194 (1987) 283. 

[25] V. L. Chernyak and A. R. Zhitnitsky, Phys. Rept. 112 (1984) 173. 

[26] V. L. Chernyak and A. R. Zhitnitsky, JETP Lett. 25 (1977) 510. 

[27] G. P. Lepage and S. J. Brodsky, Phys. Lett. B 87 (1979) 359; Phys. Rev. D 22 (1980) 2157; 
A. V. Efremov and A. V. Radyushkin, Theor. Math. Phys. 42 (1980) 97; Phys. Lett. B 94 
(1980) 245. 

[28] A. Schmedding and O. Yakovlev, Phys. Rev. D62 (2000) 116002. 
[29] I. V. Anikin, A. E. Dorokhov, L. Tomio, Phys. Lett. B475 (2000) 361. 



15 



